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AN  INTRODUCTION  TO  THE  AERODYNAMICS  OF  FLIGHT  DYNAMICS 


I 


H.  H.  B.  M./Thomas 


SUMMARY 


The  solution  of  problems  of  flight  dynamics  requires  the  aerodynamic 
forces,  which  are  called  into  play,  to  be  expressed  in  a suitable  form.  In  this 
context  a suitable  form  is  one  which  adequately  reflects  the  nature  of  the  motion 
being  considered  and  is,  at  the  same  time,  convenient  for  the  solution  of  the 
equations  of  motion.  In  the  opening  sections  of  this  paper  formulation  in  terms 
of  aerodynamic  derivatives,  and  generalizations  thereof,  are  considered.  There 
follows  a brief  discussion  in  broad  and  simple  physical  terms  of  how  Che  various 
motion  variables  give  rise  to  forces  and  moments,  which  within  a linearized 
framework  are  expressible  as  force  or  moment  derivatives,  specifically  for  an 
aeroplane,  a 


Paper  to  be  presented  as  a contribution  to  a course  of  lectures  on 
"Aerodynamic  inputs  for  problems  in  aircraft  dynamics"  at  the  von  Karman 
Institute  for  Fluid  Dynamics^  25-29  April  1977. 
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1 INTRODUCTION 

In  setting  up  a mathematical  model  for  the  calculation  of  the  dynamics  of 
an  aircraft,  it  is  necessary  to  relate  the  aerodynamic  forces  to  the  state  vari- 
ables, that  is,  those  motion  and  other  variables  necessary  to  specify  the  flight 
condition.  For  conditions  of  steady  flight,  or  conditions  of  equilibrium,  the 
problem  presents  no  inherent  difficulty  since  these  conditions  can  be  assumed  to 
have  persisted  for  a long  time.  This  renders  the  forces  fully  determinate  in 
terms  of  the  reduced  set  of  variables,  which  specify  the  steady  state.  Unsteady 
motion,  on  the  other  hand,  presents  real  difficulties  and  there  is  no  argument 
by  which,  a priori,  we  can  circumvent  the  fact  that,  strictly  speaking,  the 
aerodynamic  forces  and  moments  acting  on  an  aircraft  during  unsteady  motion  are 
'functionals',  that  is,  they  depend  on  the  whole  set  of  values  of  the  state 
parameters,  past  and  present. 

To  a large  extent  the  above  difficulty  has  been  ignored  in  the  past 
development  of  the  subject  of  flight  dynamics  and  the  formulation  of  the  aero- 
dynamic forces  and  moments  has  proceeded  along  essentially  intuitive  lines  of 

1-3 

reasoning.  Early  this  century  Bryan  introduced  a formulation  in  terms  of 
linearized  aerodynamic  derivatives,  which  expressed  the  aerodynamic  force  or 
moment  as  the  sum  of  the  steady-state  value  and  linear  terms  in  the  instantaneous 
values  of  the  disturbance  (or  perturbation)  velocities  and  motivator  angles. 

Some  terms  in  the  first  derivative  of  some  of  the  quantities  were  added  as  time 
went  on.  This  process  replaces  the  'functional'  by  a function,  in  other  words, 
if  A be  a typical  aerodynamic  reaction,  then  in  free  air  (away  from  the  Earth 
or  other  large  object)  we  may  express  the  process  mathematically  as  follows,  see 
Fig  1, 

A = A[u,v,w,p,q,r,5>n,c],  the  functional 

as  A(u,v,w,p,q,r,C»h,C,v,w,C,Tl,0 

w A + Au'+Av'  + 

e u V 


,1 

H 


where  the  suffix  e refers  to  equilibrium  conditions  of  steady  flight,  for 
example,  with  u^  and  w^^O,  v^^O,  p“q“r*0,  C“C“0  but  h “ • 

The  justi.  'on  for  the  continued  use  of  these  'quasi-steady'  aerodynamic 
derivativ  'f  the  development  of  more  generalized  formulations  stemming 

from  the  ori^x. linearized  form  is  to  be  found  in  the  fact  that  seemingly  they 
have  described  adequately  aircraft  behaviour  in  most  conditions  of  practical 


4 


Ae 

1710 


r 


interest.  Furthermore  they  have  received  a certain  degree  of  validation 
through  wind-tunnel  testing.  Nonetheless  it  must  be  emphasised  that  certain 
aspects  of  the  concept  of  derivatives  lacks  rigour.  More  recently,  Tobak^  has 
tried  to  introduce  that  rigour  by  starting  with  the  proposition  that  the  force 
or  moment  is  a functional.  He  derives  an  integral  form  for  an  aerodynamic 
coefficient  i where  k may  be  X,  Y,  Z or  £,  m,  n , which  contains 

indicial  responses  in  a functional  form.  The  subsequent  development  of  his 
formulation  hinges  on  "replacing  the  functionals  by  appropriate  functions, 
whose  dependence  on  the  past  is  denoted  by  a limited  number  of  parameters 
rather  than  by  continuous  functions".  We  cannot  go  into  the  details  of  his 
analysis  here,  but  we  note  that  he  emerges  with  formulations  akin  to  those 
discussed  later.  However  a point  worth  making  is  that  the  analysis  is 
restricted  to  motions  reproducible  in  a wind  tunnel,  namely,  a rectilinear 
flight  path  with  constant  resultant  velocity. 

Even  if  we  accept  the  derivative  formulation  with  all  its  defects  there 
remain  fairly  formidable  aerodynamic  problems  in  the  calculation  of  the  various 
quantities,  so  defined.  These  will  be  discussed  in  their  context  in  later 
sections  of  this  paper.  Furthermore  there  is,  and  as  far  as  can  be  foreseen 
always  will  be,  a need  to  have  available  methods  suited  to  the  various  stages 
of  the  design  procedure.  There  is  a willingness  to  sacrifice  some  accuracy  in 
the  interest  of  ease  of  application  and  generality  for  those  methods  to  be 
used  during  the  early  stages  of  design.  However,  the  main  aim  of  this  intro- 
ductory paper  is  to  give  in  broad  qualitative  terms  an  account  of  the  source  of 
the  various  contributions  to  the  forces  and  moments  acting  on  an  aeroplane  in 
unsteady  flight.  We  shall  not  be  concerned  with  methods  of  estimation,  per  se. 

In  order  to  achieve  the  maximum  convenience  and  generality  it  is  customary 
to  resort  to  a non-dimensional  form  for  all  the  aerodynamic  quantities  involved. 
In  this  way  the  major  effects  of  speed,  size  and  air  density  are  accounted  for 
automatically.  It  is  unfortunate  that  schemes  for  forming  these  non-dimensional 
quantities  have  not  evolved  in  an  organized  and  logic  fashion.  The  result  is 
that  the  newcomer  to  the  subject  is  faced  with  inconsistencies  which  he  finds 
difficult  to  comprehend  and  resolve.  In  what  follows, the  relationship  between 
the  coefficients  and  derivatives  in  the  system  most  commonly  used  in  the  United 
States  and  that  embodied  in  the  documents  issued  by  the  International  Standards 
Organisation  are  explained  and  detailed.  Additionally  the  corresponding 
relationship  between  the  latter  and  the  earlier  British  scheme  of  notation  is 
also  considered. 


The  use  of  the  coefficient  form  is  essentially  intended  for  treating  the 
forces  and  moments  generated  by  the  relative  motion  of  the  aircraft  and  the 
ambient  air  and  for  this  reason  the  divisors  contain  the  kinetic  pressure,  ^pV  , 
as  a factor.  For  other  forces,  even  aerodynamic  ones  such  as  the  thrust  of  an 
engine  such  divisors  are  not  entirely  satisfactory  but,  since  in  the  dynamic 
studies  it  is  necessary  to  be  consistent,  are  used.  There  is,  however,  a draw- 
back in  this  practice  as  it  can  be  misleading  in  the  process  of  forming  the  aero 
dynamic  derivatives  in  their  non-dimensional  (or  normalized)  form.  We  have  to 
draw  a clear  distinction  between  what  might  be  termed  coefficient  derivatives 
and  'coefficients  of  derivatives',  which  is  one  way  of  describing  the  usual  non- 
dimensional  aerodynamic  derivatives. 

2 CONCEPT  OF  AN  AERODYNAMIC  DERIVATIVE  AND  ITS  USES 

We  have  already  mentioned  in  broad  terms  the  concept  of  a derivative  of  a 
force  or  moment.  The  reduction  of  this  to  its  normalized  (or  non-dimensional) 
form  also  needs  consideration.  Let  us  return  to  the  equation  previously  given 
for  a typical  aerodynamic  reaction  in  the  special  case  of  motion  well  away  from 
another  body  that  can  introduce  an  interference  effect,  such  as  the  Earth.  We 
can,  without  loss  of  generality,  confine  our  attention  to  a motion  in  the  longi- 
tudinal plane  alone,  in  which  case  we  are  concerned  with  two  forces  (X,Z)  and 
one  moment  (M).  These  are  expressible,  within  the  assumption  of  small  perturba- 
tion, in  the  forms, 

X ■ Xru,w,q,ri'l  » X (u  ,w  ,ti  ) + X u'  + X w'  + X q'  + X n'  , 

L ^ J ee  e e u w q h 

Z * Zru,w,q,ril  Z (u  ,w  ,n  ) + Z u'  + Z w'  + Z q'  + Z n'  + Z.w'  , 
i-  ’ ■*  eeee  u w q n w 

and 

M = Mru,w,q,n3  « M (u  ,w  ,n  ) + M u'  + M w'  + M q'  + M n'  + M.w' 

*■  -*  eeee  u w q n w 

In  these  equations  the  quantities  are  dimensional,  although  it  is  significant 
that  a merit  of  the  new  ISO  notation  (see  also  Ref  4)  is  that  the  corresponding 
equations  in  normalized  form  have  the  same  form  exactly.  Hereafter  we  shall 
emphasise  the  difference  by  either  adding  the  superscript  0 or  the  subscript 
'ord',  when  we  wish  to  denote  a dimensional  quantity.  The  plain  symbols  are 
then  the  aero-normalized  derivatives  of  the  ISO  system,  for  example,  X^,X^,u',w 

As  a basis  for  this  normalizing  process  we  choose  divisors  formed  in  the 
same  consistent  manner  as  a system  of  units  derived  from  three  fundamental 
devisors , 


(1)  for  lengths,  (a  length  defined  by  a significant  geometric 

feature) 

(2)  for  velocities,  the  resultant  velocity  in  the  datum  flight  condition, 
V (resultant  of  u and  W in  above  example)  and 


(3)  for  forces 


ord 


From  this  choice  of  divisors  it  follows,  for  example,  that 


and 


V 


' ^ord 


ord 


with,  of  course. 


'■  ■ m. 


ord 


This  example  illustrates  the  meaning  of  the  expression  'coefficient  of  a 
derivative',  since  this  is  precisely  what  the  first  of  the  bracketed  terms  for 
Z q'  is.  It  is  perhaps  helpful  to  set  out  in  full  the  divisors  most  commonly 
needed  and  this  is  done  in  Table  1. 

It  is  implicit  in  the  form  of  the  expressions  given  above  for  the  forces 
and  moments  that  the  aerodynamic  derivative  formulation  is  applicable  only  to 
the  calculation  of  small  perturbation  motion  of  an  aircraft.  However,  it  may 
be  noted  here  that  this  refers  more  specifically  to  changes  in  the  variables 
involved  in  the  aerodynamic  terms  rather  than  others  occurring  within  the 
dynamics.  In  other  words  the  aircraft  motions  that  may  justifiably  be  calcu- 
lated on  this  basis  can  involve  appreciable  changes  in  some  state  variables, 
for  example,  the  attitude  angles.  The  derivative  form  for  the  aerodynamic 
forces  and  moments  is  particularly  well  suited  to  the  consideration  of  the 
stability  (in  the  linearized  sense)  of  a given  steady  state. 
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It  may  be  remarked  here  that  the  derivatives  used  in  the  USA  are  also 

essentially  of  the  above  type,  but  there  is  in  this  case  an  inconsistency  in 

the  unit  of  force  used  to  reduce  a force  and  a moment  to  their  respective 

normalized  forms.  For  the  force  the  divisor  is  as  above  (^p  VS)  , but  for 

2 e e ord 

the  moment  it  is  (P  V^S)  (a  representative  length).  This  usage  can  be 

e e ora  u 

traced  back  to  the  early  days  of  aviation  when  it  was  the  practice  to  use  the 

2 ...  . . 
factor  pV  in  forming  all  coefficients  and  lift  drag  and  pitching  moment 

coefficients  had  the  following  forms, 


k 


L 


L 

2 ’ 
pV  S 


2 

pV  S 


and 


k 

m 


M 


The  forms  resulting  from  the  retention  of  this  earlier  practice  for  the  moments 
is  open  to  another  interpretation  since  the  choice  made  for  the  representative 
(or  characteristic)  length  is  c/2  for  the  longitudinal  quantities  and  b/2  , 
semi-span,  for  the  lateral  quantities. 


Thus,  for  example,  we  may  write 


when  the  second  variant  of  the  expression  can  be  interpreted  as  using  c for 

4-  , in  the  ISO  notation.  This  stratagem  merely  shifts  the  inconsistency  to 

/q4o\ 

the  normalized  variables,  such  as,  — 1 and  the  time  unit, 

\ ’’e  / ord 

with  4q  “ c/2  would  not  conform  with  4^=3.  Nevertheless  the  two  sets  of 
derivatives  are  easily  related  one  to  the  other.  To  illustrate  this  we  have  to 
digress  a little  and  look  at  the  well-known  present-day  versions  of  the  aero- 
dynamic coefficients.  These  are  defined  as  follows: 


for  the  longitudinal  forces  and  moment  and 


for  the  lateral  force  and  moments,  almost  universally  and  in  particular  in  the 
USA.  We  should  remark  here  that  the  International  standard  does  not  lay  down 
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any  chosen  length  for  . However,  if  the  choice  is  made  of  adopting  * c 
for  longitudinal  quantities  and  11^=8  for  the  lateral  quantities,  we  note 
that  the  normalized  forces  and  moments  have  a form  analogous  to  the  coefficients. 
Now  the  product  of  a derivative  and  its  associated  normalized  variable  repre- 
sents the  contribution  of  that  state  variable  to  the  force  or  moment  expressed 
in  its  normalized  form.  If  we  express  this  fact  mathematically  we  have,  for 
example. 


= M w'  , by  definition 


Ae 

1710 


9 


Other  quantities  can  be  similarly  related.  Table  2 attempts  to  set  out  the 
more  or  less  complete  set  of  relationships  necessary  to  convert  from  results  in 
one  system  of  notation  to  the  other.  The  remarks  about  the  derivatives  Cx^ 
and  C2^  should  be  noted.  Also  given  in  Table  2 are  the  factors  to  be  applied 
if  the  semispan,  b/2  , is  used  in  place  of  b for  . 

3 CONCEPT  OF  A COEFFICIENT  DERIVATIVE  AND  ITS  USES 

It  not  infrequently  happens  that  a problem  in  flight  dynamics  involves 
large  perturbations  of  one  or  more  of  the  state  variables.  In  particular  there 
is  a class  of  problems  wherein  the  speed  varies  substantially.  For  such  prob- 
lems it  is  preferable  to  exclude  the  speed  from  the  expansion  of  the  aerodynamic 
forces  and  moments  by  using  what  we  shall  term  coefficient  derivatives,  that  is, 
partial  derivatives  of  the  appropriate  aerodynamic  coefficient  with  respect  to 
a variable.  Thus  typically  we  have. 


X = ^pV^SC^  : C^ 

M = ^pV^ScC  : C 
m m 


e 


C 

me 


+ C a' 
Xa 


+ C a’ 
ma 


+ C B'  + C q' 
XS  Xq^ 


+ C + C q' 
mB  mq 


+ . . . 


+ . . . 


where  the  multipliers  of  C„  and  C are  dimensional  quantities,  p being 

X m 

the  current  value  of  the  air  density  and  V the  current  value  of  the  resultant 
velocity.  In  the  expansions  of  the  coefficients  the  primed  quantities  are  aero- 
normalized  perturbation  variables.  Thus,  see  Fig  2, 


etc . , 


if  suffix  e refers  to  a synnnetrical  flight  condition.  It  should  be  noted 

that  in  this  treatment  of  the  forces  and  moments  the  current  value  of  the 

resultant  velocity,  V , is  used  in  the  definitions  and  that  all  the  coefficient 

derivatives  are  functions  of  the  current  Mach  number,  M = V/a  . This  is  in 

contrast  to  those  of  section  2 which  are  functions  of  the  datum  flight  Mach 

number.  M = V /a  . On  this  account  the  coefficient  derivatives  or  derivatives 
e e e 

of  coefficients  should  not  be  confused  with  the  usual  aerodynamic  derivatives 
used  in  the  USA,  which  as  we  have  discussed  are  essentially  of  a different  type 
in  spite  of  the  similarity  in  notation. 


It  is  possible  to  deduce  the  relationships  between  the  derivatives  of 

4 

section  2 and  those  of  the  present  section  . For  a general  datum  flight  condi- 
tion these  relationships  are  complicated,  but  they  simplify  considerably  when 
the  datum  flight  condition  takes  the  more  usual  form  of  steady,  S3rmmetric  recti 
linear  flight.  Then  for  an  arbitrary  set  of  body  axes  we  have,  for  example. 


Z 

u 


Z 

w 


cos  a - (C„  ) 
e Za  e 


sin  a + (C„  ) 
e Za  e 


sin  a 

e 


cos  a 

e 


) 


Z 

q 


(c^  ) 

Zq  e 


Z. 

w 


cos 


a 

e 


if 


2^=0  (say). 


Z 

n 


(c_  ) 

Zn  e 


where  the  suffix  e indicates  that  the  quantities  are  evaluated  at  the  datum 
condition. 


Similar  relationships  exist  for  the  moment  derivatives  except  that  the 
right-hand  sides  are  multiplied  by  an  additional  factor,  > which  can  be 

set  equal  to  unity  if  Hj  = = c (say).  For  aerodynamic  body  axes,  a^  = 0 

(Fig  3)  and  the  relationships  simplify  even  further.  In  addition  in  that  case 
we  can  introduce  the  lift  and  drag  coefficients  into  the  expressions  on  the 
right-hand  sides  (cf  section  5.1).  If  a dependence  of  upon  the  accelera- 

tion V had  been  assumed,  the  expression  for  Z^  would  have  contained  a term 
in  C • , but  it  is  normal  practice  to  ignore  this  term  just  as  we  have 

o V 

previously  ignored  terms  involving  derivatives  with  respect  to  u in  the  small 
perturbation  expressions  for  forces  and  moments  given  in  section  2. 


p 
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4 ALTERNATIVE  FORMULATIONS  OF  AERODYNAMIC  FORCES  AND  MOMENTS  FOR  OTHER 

MOTIONS 

We  have  already  remarked  in  the  introduction  that  further  generalizations 
are  needed  and  according  to  Tobak^  are  justified.  One  such  generalization  is 
that  in  which,  in  addition  to  substantial  changes  in  the  velocity,  the  changes 
in  the  angles  of  attack  and  sideslip  are  sufficiently  large  as  to  render  the 
previous  formulations  of  the  aerodynamic  terms  unacceptable.  We  then  incorporate 
within  the  first  term  of  the  expressions  for  a force  or  a moment  coefficient  the 
effect  of  the  above  two  angles  as  a function.  The  combination  of  instantaneous 
values  of  the  velocity  and  angles  of  incidence  define  a possible  steady  state, 
so  that  we  may  regard  the  expressions  for  the  aerodynamic  force  or  moment  as 
representing  a disturbance  about  the  said  steady  state. 

In  this  way  we  arrive  at  the  formulation  illustrated  by  the  following, 

and 

C = C (a, 6)  + C (a,B)n'  + C (a,6)q'  + ... 
m m mri  mq 

By  virtue  of  what  we  have  already  stated  the  first  terms  in  such  formulations 
are  obtained  from  so-called  'static'  tests  or  calculations.  Just  as  it  is 
sometimes  necessary  to  elaborate  this  formulation  further  by  including  the 
effect  of  the  motivators  (the  effect  of  q in  above  examples)  likewise  as  a 
function,  it  may  also  be  permissible  on  occasion  to  omit  the  sideslip  angle 
from  all  but  the  first  term.  The  theoretical  arguments  advanced  by  Tobak  can 
take  us  no  further  than  saying  that  strictly  the  derivatives  must  be  determined 
with  reference  to  the  state  defined  by  both  angles  of  incidence.  On  the  other 
hand  too  few  experiments  have  been  made  to  allow  of  a more  general  statement  on 
empirical  grounds. 

In  dynamic  problems  which  involve  high  rates  of  rotation,  such  as  departure 
conditions,  spin  entry  and  established  spin  motions,  further  modification  of  the 
formulation  of  the  aerodynamic  terms  is  necessary.  Again  Tobak  has  put  forward 
arguments  in  support  of  a specific  form  for  this  case,  but  we  may  approach  the 
question  from  a slightly  different  angle.  If  we  wish  to  introduce  an  angular 
velocity  into  the  first  term  of  a more  general  expression  than  those  quoted 
above,  this  angular  velocity  must  clearly  be  about  the  instantaneous  direction 
of  the  velocity,  or  along  the  tangent  to  the  flight  path,  see  Fig  2. 
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Let  us  denote  the  angular  velocity  by  u and  take  as  our  examples  two 
coefficients  which  will  be  strongly  affected  by  the  presence  of  the  angular 
velocity.  Two  such  coefficients  are  those  for  the  side-force  and  the  rolling 
moment.  These  are  now  written, 

Cy  = 0^(0, 6,0))  + C^(.0L,&,O  + C.^pP'  + CYj.r'  + ... 

and 

S = C^(a,6,o))  + C^(a,6,C)  + C^(a,6,C)  + C^^pP'  + Cj^^r'  + ... 


By  choosing  o)  such  that  it  has  a component  p , equal  to  the  current  rate  of 
rotation  about  the  x-axis,  along  that  axis  we  obtain  what  is,  in  effect,  Tobak's 
proposed  formulation^.  This  would  imply  that  the  yaw-rate  derivatives  are  to  be 
determined  (experimentally  or  theoretically)  from  motions  involving  oscillations 
about  the  coning  motion  defined  by  a,  6 and  lo  = p sec  o , see  Fig  2.  It  is, 
however,  necessary  to  note  that  the  above  choice  for  u)  is  not  the  only  one 
possible.  In  any  case  it  is  desirable  that  any  formulation  of  the  aerodynamic 
terms  be  validated  by  tests  involving  motion  of  the  most  general  kind.  It  is 
to  be  hoped  that  further  light  may  be  cast  on  the  matter  by  an  investigation, 
presently  being  undertaken,  into  the  motion  of  free-flight  models  during  post- 
stall gyrations,  spin  entry  conditions  and  in  established  spins  of  both  the 
steep  and  flat  type. 


5 ON  HOW  THE  FORCES  AND  MOMENTS  ARISE 


In  this  section  the  aim  is  to  describe  in  the  broadest,  and  hopefully 
simplest  possible,  terms  the  nature  of  each  contribution  to  the  overall  forces 
and  moments  and  to  identify  those  components  of  an  aeroplane  which  make  signifi- 
cant contributions  to  the  individual  derivatives. 


It  is  desirable  to  work  in  some  specified  system  of  axes  and  the  aero- 
dynamic-body axes  are  chosen  for  their  convenience  and  general  nature.  A 
definition  of  this  axis  system  is  not  amiss,  so  we  note  that  aerodynamic-body 
axes  are  such  that  in  the  datum  flight  condition  the  x-axis  is  parallel  to  the 
projection  on  to  the  plane  of  symmetry  of  the  velocity  of  the  origin  of  the 
axes  (usually  the  centre  of  gravity  of  the  aeroplane).  The  x-axis  is  positive 
towards  the  nose  of  the  aeroplane,  the  z-axis  is  parallel  to  the  plane  of 
symmetry,  positive  ventrally.  The  y-axis  is  positive  to  starboard;  when  the 
centre  of  gravity  lies  in  the  plane  of  symmetry,  the  latter  is  the  zx-plane, 
see  Fig  3. 
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The  choice  of  axes  made  here  does  not  imply  that  these  axes  are  the  most 
suitable  for  all  aircraft  configurations  or  for  all  experimental  techniques. 

It,  therefore,  sometimes  becomes  necessary  to  be  able  to  convert  the  \alues  of 
the  aerodynamic  derivatives  with  respect  to  one  system  of  axes  to  those  appropri- 
ate to  another  system  of  axes.  This  matter  is  discussed  in  section  6. 

To  avoid  undue  complexity  we  have  previously  excluded  from  the  discussion 
any  effect  location  of  the  aircraft  (that  is,  the  position  of  its  centre  of 
gravity  and  its  attitude  angles  or  orientation)  may  have  on  the  aerodynamic 
forces  acting  upon  it.  This  is  not  always  permissible.  For  example,  the 
proximity  of  an  aeroplane  to  the  ground  (Earth's  surface)  as  during  landing  and 
take-off  introduces  changes.  Six  parameters  govern  these  effects,  namely,  x^, 
y^,  ZqJ  f,  0,  <P,  of  which  the  first  three  define  the  position  of  the  aeroplane's 
centre  of  gravity  with  respect  to  an  Earth-fixed  system  of  axes  and  the  second 
trio  the  aeroplane's  orientation  relative  to  these  axes.  Since  it  is  usually 
sufficient  to  simplify  to  a flat  Earth  the  first  trio  may  be  replaced  by  h , 
the  height  of  the  centre  of  gravity  above  the  ground.  Wind-tunnel  test  results 
as  well  as  those  from  theory  have  indicated  that  some  of  the  derivatives  already 
discussed  are  subject  to  considerable  modification.  However,  some  care  needs  to 
be  exercised  in  interpreting  results  of  this  sort  as  they  not  infrequently  refer 
to  a constrained  flight  condition.  In  the  usual  wind-tunnel  test  the  centre  of 
gravity  remains  fixed,  that  is,  no  change  in  flight  path  occurs.  Orientation 
of  the  aeroplane  with  respect  to  the  incident  flow  becomes  indistinguishable 
from  orientation  with  respect  to  the  Earth.  To  illustrate  this  point  it  is 
convenient  to  refer  the  motion  to  a suitable  body  axis  system,  see  Fig  3,  for 
example,  with  the  x-axis  aligned  with  the  no- lift  line  of  the  aeroplane.  In 
Fig  4 two  special  motions  of  the  aeroplane  are  shown,  one  in  which  the  attitude 
in  pitch  is  zero  (or  constant)  and  the  angle  of  attack  varies  the  other  in  which 
the  angle  of  attack  is  zero  whilst  the  attitude  changes.  These  and  the  more 
general  motion  shown  in  Fig  5a  may  be  contrasted  with  the  flight  conditions 
usually  examined  theoretically  and  experimentally.  Fig  5b.  Whether  or  not  these 
are  significant  derivatives  with  respect  to  the  attitude  angle  © has  yet  to  be 
determined.  Similar  remarks  apply  to  the  other  two  attitude  angles,  'i  and 

The  same  diagrams  serve  to  illustrate  the  essential  difference  between  the 
aerodynamic  forces  arising  from  the  rate  of  pitch,  q and  those  due  to  the  rate 
of  change  of  the  angle  of  attack,  a . Again  there  is  parallel  to  be  drawn  for 
the  rate  of  yaw,  r and  the  rate  of  change  of  sideslip  angle,  8 . 


5.1  The  longitudinal  X-force  derivatives,  X and  X 


w 


The  derivative  X^  represents  the  contribution  to  the  force  along  the 

x-axis  due  to  unit  incremental  velocity  along  the  x-axis,  so  that 

incremental  force  due  to  u'  . Likewise  X w'  is  the  corresponding  increment 

w 

in  force  due  to  w'  . Since  lift  and  drag  coefficient  data  are  more  familiar 

and  generally  more  readily  available,  we  shall  express  these  derivatives  X^ 

and  X in  terms  of  derivatives  of  these  coefficients, 
w 

2 

We  may  write  X = ^pV  SCj^(M,R,a,n,q,d)  as  the  dimensional  form  of  the 
X-force,  from  which  we  have  the  derivative 


where 


Y represents  each  of  the  parameters  within  the  brackets  for  the  function 
C.^  and  the  subscript  e denotes  that  the  whole  is  evaluated  at  datum 
conditions.  Now  in  general 

„2  2^2  2 

V = u + V + w , 


tan  a = — 


M = - 


R = 


pay 

V 


where  we  have  for  ease  of  writing  omitted  the  superscript  0 or  suffix  ord  . 
For  small  perturbations  in  u,  v,  w and  in  a these  can  be  written,  for  the 
chosen  axes. 
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in  the  datum  flight  conditions.  Evaluation  of  the  derivative  in  its  normalized 
form  now  gives 


By  virtue  of  the  fact  that  to  first  order  in  a'  we  may  write  C = - Co'  - C , 

A Li  U 

the  above  expression  can  be  expressed  in  terms  of  the  drag  coefficient,  thus, 


The  third  term,  which  represents  the  effect  of  the  Reynolds  number,  is  often 
ignored,  since  it  is  assumed  that  either  the  estimation  or  measurement  of  the 
drag  has  been  made  at  a representative  scale.  In  deriving  this  result  for 
we  have  restricted  our  attention  to  airframe  aerodynamic  forces,  which 
necessitates  a separate  treatment  of  the  thrust  contribution  to  the  X-force. 

This  course  is  advisable  since  the  introduction  of  a thrust  coefficient  can  be 
inconvenient  and  even  misleading.  If  T represent  the  thrust  and  i^  the 

inclination  of  the  thrust  line  to  the  x-axis,  we  have  components  of  thrust  along 

T T . . 

the  x-axis  and  z-axis  of  X = T cos  and  Z = T sin  respectively.  The 

contribution  of  the  thrust  to  X can  be  written 

u 


u 


It  is  clearly  zero  in  gliding  flight  since  then  T * 0 , otherwise  we  require 
to  know  the  dependence  of  the  engine  thrust  on  the  aeroplane  speed.  For  example, 
if  the  horsepower  of  the  engine  is  constant,  TV  = constant,  so  that 

/0T  \ 'j' 

I- — r)  “ ■ TT"  » which  is  a reasonable  approximation  in  the  case  of  a piston- 

you  /g  Vq 

engined  aeroplane  with  variable-pitch  propellers.  There  is  much  less  dependence 
of  thrust  on  speed  for  jet  or  rocket  engines . 


X 

u 


U U 


In  general 


where  the  Reynolds  number  term  has  been  ignored.  By  application  of  the  same 
basic  relationship  of  force  and  coefficient  and  evaluating  the  derivative  in  a 
similar  way  we  have 

i) 

w Xa ' e 

which  by  virtue  of  the  relationship  between  C^,  and  can  be  rewritten 


X = X 
w w 


Ct  - ,) 

L Da'  e 

e 


Estimates  of  the  derivatives  involved  in  the  above  expressions  for  X^  and  X 
can  usually  be  obtained  from  generalized  data  and  are  also  available  from  wind 
tunnel  test  data  obtained  at  an  early  stage  in  the  design  procedure. 

5.2  The  Z-force  derivatives.  Z . Z and  Z 


The  following  results  are  readily  obtained  by  applying  arguments  on  the 


lines  of  those  of  section  5.1. 


"u  V 

e e 


S-.&).  ^ 

S \3u/ 


which  when  combined  yield 


Z = Z - 20- 
u u L 


"e  (&l 


of  which  the  first  term  is  usually  so  small  that  it  is  ignored. 
For  the  derivative  with  respect  to  w'  we  have, 


3a  / ’ 

' e e 


where  this  time  the  second  term  is  usually  small  compared  with  the  first  term. 

As  for  the  X force  derivatives  the  information  on  the  lift  and  drag 
coefficients  is  so  basic  to  the  design  of  the  aeroplane  that  the  evaluation  of 
the  Z force  derivatives  also  presents  no  problems. 
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The  derivative  Z is  associated  with  the  angular  velocity  of  the  aircraft 

q 

in  pitch.  We  have  already  seen  that  motions  can  be  imagined  in  which  there  is 
rotation  in  pitch  without  change  in  the  angle  of  attack  as  well  as  motion  in 
which  the  angle  of  attack  varies  without  any  rotation.  Accordingly  q and  d 
are  variables  with  their  separate  effects  just  as  3 and  a are.  Reference 
to  Figs  6 and  7 shows  that  the  effect  of  the  rotation  rate  q about  the  y*axis 
is  equivalent  to  a decrease  in  the  effective  angle  of  attack  of  those  portions 
of  the  aircraft's  body  and  wing  lying  forward  of  the  y-axis  and  an  increase  in 
the  effective  local  angle  of  attack  of  the  tail  and  those  portions  of  the  wing 
and  body  behind  the  y-axis.  Within  the  assumptions  of  linearized  theory  there 

is  a strict  equivalence  in  the  load  distribution  of  the  cambered  wing  in  recti- 

linear flight  and  the  rotating  wing.  The  pressure  distribution  is  proportional 

to  q/V  for  the  contributions  from  the  wing,  body  and  tailplane.  The  integral 
e 

of  these  pressure  distributions  gives  the  overall  force  and  hence  its  derivative 
. A simple,  but  approximate,  expression  for  the  tailplane  contribution  can 
be  obtained  by  ignoring  the  effect  of  the  additional  downwash  induced  by  the 
additional  lift  distribution  on  the  wing  due  to  the  rate  of  rotation,  q . We 

further  replace  the  distribution  of  local  effective  incremental  angle  of  attack 

q'^T  ... 

over  the  tailplane  by  a mean  value,  — . Here  is  the  tail-arm  and  is 

e 

usually  taken  as  equal  to  the  distance  from  the  centre  of  gravity  to  the  aero- 
dynamic centre  of  the  tailplane,  see  Fig  7. 

The  mean  kinetic  pressure  over  the  tailplane  is  affected  by  the  wing  wake 
associated  with  the  combined  a and  q loadings.  Let  the  ratio  of  this  mean 
to  free-stream  kinetic  pressure  be  , then  we  have  approximately. 


(Z  ),  , = - V Q S 

t due  to  q e e^q 


which  yields 


^ ‘it  • 

The  distribution  of  effective  angle  of  attack  due  to  the  pitch  rate  also 
produces  a contribution  to  the  rotary  derivative  Z^  from  the  body  and  the 
wing.  That  of  the  body  can  be  estimated  using  slender-body  theory  or  one  of 
the  more  refined  theoretical  methods.  Strictly  speaking  allowance  should  be 
made  for  the  mutual  interference  of  wing  and  body,  but  this  is  seldom  done, 
because  of  the  more  dominating  effect  of  the  tailplane.  Likewise  the  cancelling 
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effects  of  reduction  of  effective  angle  of  attack  of  that  part  of  the  wing 
forward  of  the  centre  of  gravity  and  increase  in  effective  angle  of  attack  of 
those  parts  aft  of  the  centre  of  gravity  tend  to  make  the  wing  contribution 
small,  except  for  highly-swept  wings.  However,  there  is  no  inherent  difficulty 
in  applying  any  of  the  available  lifting  surface  theories  to  the  calculation  of 
the  derivative  Z for  a wing  alone.  It  is  usually  accepted  that  simple 
summation  is  adequate  for  yielding  the  value  of  the  derivative  for  the  complete 
aircraft.  This  view  is  further  justified  by  the  fact  that  this  particular 
derivative  is  of  little  importance  in  the  dynamics  of  most  aeroplanes,  apart 
perhaps  for  those  having  very  small  values  of  the  relative  density  parameter, 

h • 

5.3  Moment  due  to  change  in  velocity  component  along  x-axis  - 

Just  as  we  expressed  the  force  derivatives  and  Z^  in  terms  of  the 

more  familiar  lift  and  drag  coefficients  so  we  may  express  the  moment  deriva- 
tive, M , in  terms  of  a derivative  of  C , at  any  rate  as  far  as  the  contribu- 
u m 

tion  from  the  airframe  forces  is  concerned.  Accordingly  we  write 


= ^pV^ScC  (M,R,a,n,q,d) 


as  an  equation  for  the  dimensional  moment  in  terms  of  its  coefficient. 

From  the  above  relationship  we  again  obtain  on  differentiation  the 
dimensional  derivative 

where  as  before  y represents  each  of  the  parameters  within  the  bracket  for 

the  function  C and  the  subscripts  'ord*  and  e denote  dimensional  value  and 
m 

datum  conditions  respectively.  Evaluation  of  this  expression  using  the 
relationships  given  in  section  5.1  yields 


/SC  SC  \ 


If  the  thrust  line  passes  above  the  centre  of  gravity  there  will  be  another 
contribution  to  . Since  the  moment  arm,  d , is  usually  small  in  terms  of 
the  mean  chord  the  contribution  is  usually  small.  We  have 
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and  so 


In  the  absence  of  compressibility  and  Reynolds  number  effects,  ■ 0 for  a 
rigid  aircraft  in  gliding  flight. 

5.4  Moment  derivative  due  to  the  change  in  velocity  component  along  the 

z-axis  - M 
w 

We  can  once  more  approach  the  estimation  of  the  moment  derivative  due  to 
the  velocity  component  perturbation,  w'  via  a relationship  between  it  and 
a derivative  of  the  moment  coefficient.  It  is  unnecessary  to  repeat  the  argu- 
ment in  detail  {af  5.2)  and  we  content  ourselves  with  merely  quoting  the  result 


In  an  analysis  of  the  pitching  moment  produced  by  changes  in  the  angle  of  attack 
it  is  customary  to  reduce  the  system  of  forces  acting  on  the  wing/body  combina- 
tion to  a moment  about  an  axis  through  the  aerodynamic  centre  and  forces 

Lq  (lift)  and  (drag)  passing  through  this  point.  For  simplicity  we  shall 

ignore  the  small  terms  arising  from  the  drag  force  and  the  thrust.  It  then 
follows  that  if  is  the  lift  on  the  tailplane  the  pitching  moment  about  an 

axis  through  the  centre  of  gravity  can  be  written, 

M - M-  + L-(h  - h_)c  - - (h  - h-)c| 

0 0 0 t ' a.c  O’ 

where  he  denotes  the  position  of  centre  of  gravity  on  the  mean  aerodynamic 
chord, 

h^c  the  position  of  the  aerodynamic  centre, 

and  I the  distance  between  the  aerodynamic  centre  of  wing-body  combination 
and  that  of  the  tailplane. 

But  the  total  lift  L * , so  that 

M = M-  + L(h  - h-)c  - LI 

0 0 t a.c 


r 
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Now 


L * pV'^S  a o 


where  a 


t eff 


“ (i  • If)  * \ • 


Hence 


C = C 
m 


/S  «,  \ 

+ (h  - h-)C  - Q a a 

“O  0 ^ \ Sc  / 


eff 


and 


m.  ■ ■■  - V (&).  ■ .„..  e^)  (.  - a. . 

Here  we  have  made  the  usual  assumption  that  the  aerodynamic  centre  position  is 
insensitive  to  variation  in  the  angle  of  attack  in  accord  with  the  findings  of 
linearized  theory. 

The  estimation  of  the  derivative  M , which  determines  the  increment  in 

w 

pitching  moment  due  to  a perturbation  velocity,  w'  , thus  reduces  to  the 
estimation  of  the  aerodynamic  centre  of  the  wing-body  combination,  the  rate  of 
change  of  the  lift  coefficient  with  angle  of  attack  together  with  the  downwash 
at  the  tailplane  and  the  lift  coefficient  derivative  for  the  tailplane,  a^^^  . 
This  last  should  strictly  take  account  of  interference  between  the  tailplane  and 
the  rear  body  of  the  aircraft. 

5.5  Rotary  damping  in  pitch  or  the  moment  derivative,  M 

— — ■■  a 

We  have  already  demonstrated  that  the  rate  of  pitch  can  be  treated  as 
equivalent  to  an  incremental  change  in  the  angle  of  attack  of  various  parts  of 
the  airframe.  The  forces  associated  with  these  incremental  angles  of  attack 
(a  longitudinal  camber)  give  rise  to  moments  about  the  y-axis  and  hence  a 
derivative  M 

q 

In  dimensional  form  the  moment  due  to  the  tailplane  is 


t due  to  q 


^^t^due  to  q ^t 


I 
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which  by  virtue  of  the  relationship  between  the  Z-force  and  q , previously 
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given,  can  be  expressed  in  derivative  form  as 


(M  ) 

q t 


) Ml 

SS2  "it 


In  contrast  to  what  happens  in  the  case  of  the  force  derivative  the  moment  due 
to  an  elementary  area  of  the  wing  ahead  of  the  y-axis  adds  to  that  contributed 
by  an  elementary  area  aft  of  the  y-axis,  so  that  the  wing  contribution  to 
can  be  important,  particularly  so  for  highly  swept  wings.  Similar  remarks 
apply  to  the  contribution  of  the  body  and  these  contributions  can  be  estimated 
on  the  same  basis  as  outlined  in  the  discussion  of  Z 

q 

5 .6  The  w or  acceleration  derivatives 

In  the  sense  that  the  derivatives  we  have  thus  far  discussed  can  be 
identified  with  forces  and  moments  acting  during  a steady  motion  (a  given  angle 
of  attack  or  angular  velocity)  they  differ  from  those  that  arise  from  a changing 
component  of  velocity  along  the  z-axis,  or  w . To  differentiate  them  they  are 
often  termed  quasi-steady. 

Consider  what  happens  to  the  pressure  distribution  over  a wing  under 
conditions  of  varying  angle  of  attack.  It  takes  time  to  adjust  and  so  differs 
from  the  pressure  distributions  associated  with  the  sequence  of  angles  of 
attack.  The  changing  conditions  are  usually  accounted  for  by  introducing  an 
acceleration  derivative  into  the  force  and  moment  formulation. 

It  can  be  demonstrated  that  this  procedure  is  not  wholly  satisfactory, 
in  general.  In  particular  if  the  angle  of  attack  of  a three-dimensional  wing 
is  suddenly  changed  from  zero  to  a fir  /alue,  the  lift  generated  will  vary 
and,  apart  from  a very  short  interval  o..  time  in  the  case  of  subsonic  flow,  is 
below  (or  lags  behind)  the  asymptotic  value,  which  is  reached  in  a finite  time 
supersonically,  but  after  a theoretically  infinite  time  subsonically . Again  if 
we  consider  a motion  in  which  q > 0 and  d = constant  it  can  be  shown  that 
the  lift  associated  with  the  changing  angle  of  attack  is,  according  to  linear- 
ized theory,  as  shown  in  Fig  9. 

A more  usual  approach  to  the  determination  of  the  w derivatives  is  via 
the  aerodynamics  of  wings  (or  complete  aeroplanes)  in  oscillatory  flow  retaining 
only  first  order  terms  in  the  frequency  parameter.  In  experiments  with  oscillat- 
ing models  the  derivatives  M and  M.  are  obtained  in  combination  as  the  sum 

q w 

(M  M.)  , but  either  a theoretically  determined  value  of  M (or  a separate 
q w q 
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steady  test)  enables  the  derivative  to  be  isolated.  The  same  oscillatory 

theory  can  be  used,  in  principle,  to  evaluate  the  tailplane  contribution  to  the 

derivative  M.  . However,  some  insight  can  be  obtained  by  an  approximate 
w 

approach,  which  accounts  for  the  lags  in  the  arrival  at  the  tailplane  of  the 
downwash  due  to  the  loading  associated  with  the  angle  of  attack. 


Convection  of  the  vorticity  downstream  implies  that  under  changing  angle- 
of*attack  conditions  the  change  in  circulation  around  the  wing  is  felt  at  the 
tailplane  after  a time  lapse  of  instantaneous  downwash  at  the 

tail  is  that  which  corresponds  not  to  the  current  angle  of  attack,  o(t)  but 


The  associated  change  in  downwash  is  then,  approximately. 


Ae 


3e  t 
3a  V 


Aa  , the  change  in  effective  angle 
of  attack  of  the  tailplane. 


This  generates  a force  and  a pitching  moment  which  can  be  written 


Z,  - - iq^pv  , 

and 

• 

From  which  the  contributions  to  the  derivatives  are  easily  obtained.  They  are. 


(Z.)^ 
w t 


Q 

^a  It 


— (—] 


and 


(M. 
w t 


Vit 


t t / 3e  \ 

“2  130  / 
Sc  ' 'e 


5.7  Control  derivatives  - Z and  M 

n n 

Deflection  of  the  tailplane,  or  part  of  it  (elevator),  with  respect  to 
the  remainder  of  the  aircraft  produces  a change  in  circulation  around  the  tail- 
plane  and  a force.  This  force  is  the  source  of  the  derivative  Z^  and  by  the 
moment  it  generates  around  the  y-axis,  the  derivative  . We  have  already 
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This  derivative  is  analogous  to  , but  in  this  case  there  is  no  large 
aerodynamic  surface  like  the  wing  involved.  When  an  aeroplane  is  in  sideslip 
components  of  it  are  at  an  angle  of  incidence  to  the  velocity  vector  (or  to  the 
relative  incident  flow),  see  Fig  2,  It  can  be  readily  appreciated  that  the  main 
contributors  to  the  sideforce  derivative  are  the  fin  and  the  fuselage.  The 
first  of  these  is  a ' lifting ' -surface,  which  can  be  accorded  an  'effective  lift 
coefficient  derivative',  a^^^  in  the  same  way  as  we  introduce  a^^  in 
section  5.  By  following  the  same  line  of  argument  we  can  easily  obtain  the 
result,  see  Fig  12. 


where  a is  the  sidewash  induced  at  the  fin  by  the  asymmetric  loading  distribu* 
tion  present  over  the  wing-fuselage  combination.  Q is  a relative  kinetic 

p 

pressure  factor,  which  except  in  unusual  circumstances,  large  angle  of  attack 
or  fin  immersed  in  a propeller  slipstream,  does  not  differ  radically  from  unity. 

6.2  Rolling  and  yawing  moments  due  to  sideslip  - L and  N 

V V 

Three  geometric  features  of  the  wing  influence  the  rolling  moment  acting 
on  an  aircraft  in  sideslipping  motion.  These  are  its  dihedral  angle,  its  posi- 
tion on  the  body  section  (high,  mid  or  low  wing)  and  its  sweep  angle.  We 
consider  these  in  turn. 

To  examine  the  effect  of  dihedral  it  is  convenient  to  restrict  our 
attention  to  a planar  wing  (see  Fig  10),  in  which  the  port  and  starboard  portions 
of  the  wing  are  inclined  at  a small  constant  angle,  T , to  the  plane  containing 
the  root  chord  and  the  y-axis. 

As  can  be  seen  from  the  diagram.  Fig  10,  the  velocity  normal  to  the  star- 
board panel  is  w^  + vT  , whilst  that  of  the  port  panel  is  w^  - vF  . This 
antisymmetric  change  in  the  effective  angle  of  attack  produces  a load  distribu- 
tion giving  a negative  rolling  moment  proportional  to  vF  , hence  a contribution 
to  the  derivative  proportional  to  the  dihedral  angle. 

To  explain  the  way  the  wing  position  on  the  body  affects  we  use  the 

concepts  of  slender-body  theory.  This  reduces  to  a consideration  of  the  two- 
dimensional  flow  associated  with  the  cross-flow  component  of  velocity  alone, 
that  is,  the  component  v , when  = 0 . Fig  11  illustrates  the  flow  around 
the  body  and  this  shows  that  the  body  induces  velocities  which  enhance  the 
dihedral  effect  when  the  wing  is  set  high  on  the  body.  On  the  other  hand  a low 
wing  will  diminish  the  dihedral  effect. 

Any  attempt  at  calculating  the  rolling  moment  due  to  sideslip  for  a swept 
wing  with  zero  dihedral  presents  a number  of  problems.  Nevertheless,  the  nature 
of  the  rolling  moment  associated  with  sweep  can  be  appreciated  using  simple 
physical  arguments.  Sideslip  (see  Fig  12)  effectively  decreases  the  sweepback 
of  the  starboard  panel  and  increases  the  sweepback  of  the  port  resulting  in  an 
increase  of  lift  on  the  starboard  side  and  decrease  on  the  port  side  so  giving 


Ae 

1710 


25 


rise  to  a rolling  moment.  Further  insight  into  the  form  of  the  rolling  moment 
can  be  obtained  if  we  consider  a large  aspect  ratio  wing  of  constant  chord  and 
sweepback,  A , see  Fig  12. 

Suppose  represents  the  component  of  normal  to  the  leading  edge 

and  the  component  of  parallel  to  the  leading  edge  of  the  wing.  Then 

V = V cos  (A  - 6)  and  V = V sin  (A  - S) 
n e t e 

at  some  point  of  the  starboard  panel.  Only  the  former  component  produces  a 

2 

lift  and  on  an  elementary  strip,  PQ,  , this  is  V a a c 6s,  where  a is 

^ r.  ennnn  n 

the  lift  coefficient  derivative  of  the  'normal'  section  and  c , 6s  are  as 

n 

marked  in  the  diagram. 

The  velocity  normal  to  the  wing  panel  is 

w = aV  = ua  = Va 
n n e e 

where  the  angles  of  attack  are  measured  relative  to  the  no-lift  line.  Only  the 
velocity  component  produces  a lifting  force  and  the  lift  on  an  elementary 

strip  normal  to  the  wing  leading  edge  is 


6L  = hp  V a a c 6s 
e n n n n 


so  that 


= V V a ac6y  , see  Fig  12. 
e n e n 


V = V (cos  A + 3 sin  A)  to  first  order  in  3 
n e 


6L  = ^p  V a (cos  A + 3 sin  A)ac6y  , 
e e n 


= ^p  V a(l  + 3 tan  A)ac6y  , 
e e 


where  a is  lift  coefficient  derivative  for  'normal'  section  and  a is  the 
n 

corresponding  quantity  for  a streamwise  section. 


Since  a corresponding  section  on  the  port  wing  experiences  a decrease  in 
the  lift  due  to  3 , a rolling  moment  is  produced  which  is  proportional  to  -o3  . 
Hence  this  particular  contribution  to  the  rolling  moment  due  to  sideslip  cannot 
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be  expressed  as  a linear  derivative  except  when  the  angle  of  attack  remains 

unchanged,  as,  of  course,  is  the  case  within  the  linearized  lateral  motion. 

The  moment  about  the  x-axis  of  the  fin  sideforce  due  to  sideslip  makes 

another  contribution  to  the  derivative  L . If  -z^  denote  the  z-coordinate 

V fv 

of  the  aerodynamic  centre  of  the  load  on  the  fin  arising  from  the  sideslipping 
motion  and  the  sideforce,  the  rolling  moment  is  , see  Fig  12.  By 

virtue  of  the  relationships  already  quoted  for  we  then  have 


L 


vf 


®lf^v 


9a  \ ^f ^fv 
96  /g  Sb 


Ae 

1710 


The  same  sideforce  also  produces  a yawing  moment  and  this  involves  the 
fin  arm,  , in  place  of  easily  shown  that 

N - n A i£\ 

~ ®lf^v  \ 96  ) Sb 

'e 

This  positive  contribution  to  the  of  the  aircraft  is  opposed  by  an  often 

large,  negative  contribution  from  the  body.  Inviscid  flow  theory  and,  in 
particular,  slender-body  theory  accounts  for  such  a contribution,  but  as  these 
theories  do  not  give  rise  to  a sideforce  the  moment  is  also  in  error. 

6.0  Rolling  and  yawing  moment  derivatives  due  to  the  rate  of  roll  - L and  N 

P.-  . . -P 

Suppose  the  aircraft  rolls  with  rate  of  roll,  p , about  the  aerodynamic- 
body  axis.  Ox  , as  shown  in  Fig  13.  This  motion  clearly  induces  a downward 
velocity  equal  to  py  for  all  points  of  the  section  of  the  starboard  wing  at  a 
distance  y from  the  plane  of  symmetry.  In  the  presence  of  the  forward  velo- 
city this  downward  velocity  can  be  interpreted  as  equivalent  to  an  increase 

in  the  local  angle  of  attack  of  see  Fig  13.  There  is  an  equivalent 

decrease  in  the  angle  of  attack  of  the  corresponding  section  of  the  port  wing. 
We,  therefore,  can  regard  a rolling  velocity  as  equivalent  to  an  antisymmetric 
twist,  which  increase  in  magnitude  to  a maximum  at  the  wing  tips.  There  is  an 
associated  antisymmetric  loading  and  a rolling  moment  proportional  to  p . 

This  is  the  source  of  the  wing  contribution  to  L . Smaller  contributions 

P 

come  from  the  tail  surfaces.  That  due  to  the  tailplane  is  the  result  of  the 
antisymmetric  twist  equivalent  of  the  rolling  and  the  changes  this  effectively 
produces  in  the  tailplane  angle  of  attack. 


■f ' ' MP  fm 
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As  for  other  fin  contributions  the  can  be  estimated  on  the  basis  of 

a mean,  effective  change  in  the  angle  of  incidence  of  the  fin.  This  can  be 
approximated  as 


- i£ 


V3  \v 


in 

which  gives  rise  to  a sideforce,  see  Fig  13, 


9a 


IffE  . 

^e  p 


(^)J 


'if  = - V Q a_6pb  ^ 

f e e p If  1 


r 

From  which  we  can  derive  the  fin  contribution  to  the  rolling  moment  derivat 


ive 


as 


S.2-  I 

Q a -XiE  J 

> If  Sb  ^ 


-i£  - 9° 

b 


TiiA 

U I 

\ e/. 


1, 


Associated  with  the  as3mnnetric  loading  that  gives  rise  to  the  derivative  L 
are  changes  in  the  forces  induced  in  the  xy-plane  and  these  produce  a yawing 
moment  proportional  to  the  rate  of  roll  and  hence  a derivative,  N 

P 

Again  there  is  a fin  contribution  and  this  can  easily  be  shown  to  be 
given  by  the  relationship, 


pf 


Vif 


Vf 

Sb 


In  both  the  expression  for  and  the  distance  represents  the 

distance  between  the  x-axis  the  aerodynamic  centre  of  the  fin  loading  due  to  p 
that  is,  due  to  a linear  twist  increasing  towards  the  fin  tip. 

6.4  Rolling  and  yawing  moment  derivatives  due  to  rate  of  yaw,  L and  N 

• — r r 

When  an  aircraft  experiences  a rate  of  yaw  a section  of  the  wing  on  the 
starboard  side  at  a distance  y from  the  plane  of  symmetry  has  a velocity 
parallel  to  the  x-axis  of  - ry  , whilst  the  corresponding  section  of  the 
port  side  has  an  increased  velocity,  + ry  (Fig  14).  This  asymmetric 
velocity  increment  decreases  the  forces  on  the  starboard  wing  and  increases 
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those  acting  on  the  port  wing.  The  resulting  moments  are  proportional  to  the 
rate  of  yaw  and  account  for  the  wing  contributions  to  and  . 

There  is  a significant  contribution  to  from  the  fin.  A rate  of  yaw 

r produces  an  effective  increment  in  the  fin  angle  of  incidence,  see  Fig  14  of 


rb 

V 

e 


which  gives  rise  to  a sideforce  on  the  fin  of 


(Y.). 

f due  to  r 


hp  V^s  Q a,  - 
e e f r If  V 

e 


and  in  turn  yields  a derivative  contribution, 


Since  we  are  concerned  with  a mean  angle  of  incidence  z,.  z, 

fr  fv 

There  is  an  even  larger  contribution  to  , which  can  be  easily  shown 
to  be  given  approximately  by  the  relationship. 


There  is  a body  contribution  to  , which  is  analogous  to  its  contribution 
to  M . 

q 

6.5  Lateral  and  directional  control  derivatives 

The  forces  and  moments  acting  on  an  aircraft  due  to  deflection  of  its 
lateral  and  directional  motivators  depend  upon  the  form  these  motivators  take. 

To  illustrate  the  derivatives  involved  we  consider  the  common  flap  type 
motivators,  that  is,  ailerons  and  rudders. 

Deflection  of  a trailing-edge  flap  alters  the  local  circulation  around 
the  aerofoil.  This  effect  extends  the  span  of  motivator  and  beyond.  In  general 
forces  and  moments  are  produced  by  the  resulting  changes  in  loading  distribution. 
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For  the  conventional  aileron  layout  (see  Fig  15)  the  downgoing  starboard  aileron 
increases  the  lift  on  the  outer  portions  of  the  starboard  wing,  whilst  on  the 
port  side  the  upgoing  aileron  decreases  local  lift.  A negative  rolling  moment 
is  thus  produced  and  associated  with  the  induced  effects  of  the  asymmetrical 
loading  and  yawing  moment.  Both  these  moments  are  proportional  to  the  mean 
angle  according  to  linearized  theory.  The  yawing  moment  is  usually  positive 
(that  is,  in  the  adverse  yaw  sense).  The  moments  are  represented  by 
N^C'  for  an  increment  C'  in  the  aileron  angle. 

Deflection  of  the  rudder  to  port  (positive  sense)  results  in  a positive 
sideforce  and  hence  a yawing  moment  (the  primary  function  of  the  rudder)  in  the 
negative  sense.  The  rate  of  change  of  this  moment  with  the  rudder  angle,  C , 
defines  the  derivative  . Because  the  aerodynamic  centre  of  the  loading  on 
the  fin  is  generally  off  the  x-axis  a rolling  moment  is  produced  at  the  same 
time  as  the  yawing  moment. 

if  the  chord-ratio  of 
v 

the  rudder  does  not  vary  much  across  the  fin. 


We  may  take  to  be  approximately  equal  to  , 


Corresponding  to  the  fin  effective  lift  coefficient  slope,  ®j£.  there  is 


a rudder  parameter  a^^  and  so  we  may  easily  deduce  the  results 

N.  = - Q„ 


^f^f 


Sb  2f 


and 


= Q, 
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7 FORCE  AND  MOMENT  DERIVATIVES  FOR  DIFFERENT  SYSTEMS  OF  AXES 

We  have  already  remarked  that  a set  of  derivatives  is  particular  to  the 
choice  of  axis  system.  It  is  frequently  necessary  to  convert  the  values  of 
derivatives  to  correspond  to  a new  choice  of  axes.  The  most  common  requirement 
is  to  convert  from  one  body  system  to  another  body  system  having  parallel  y-axis 
Let  the  two  systems  of  axes  be  Oxyz  and  • Further  let  the  origin 

0^  have  coordinates  h,0,k  relative  to  0 , in  the  system  Oxyz  , whilst  0 
has  coordinates  -hj^,0,-k^  relative  to  0^^  in  the  system  (see 

Fig  16).  Then,  if  we  write  c “ cos  e and  s “ sin  e , where  e is  the  angle 
between  taken  as  positive  if  clockwise  rotation  about 

the  y-axis  through  this  angle  brings  Ox  to  a position  parallel  to  O^x^  . 
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The  following  relationships  then  apply, 
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ch  - sk  ; 

h 

Chi 

+ skj^ 

>'1  = 

ck  + sh  ; 

k 

B 

cki 

- shj 

velocity  components. 

Pi  “ 

cp  - sr  ; 

P 

= 

CPi 

+ srj^ 

'll  = 

q ; 

q 

S 

^1 

"l  “ 

cr  + sp  ; 

r 

= 

- sp^ 

If  we  define  u,  v,  w as  the  components  of  the  velocity  of  0 along  the  axes 
X,  y and  z whilst  v^,  w^^  are  the  components  in  the  Oj^Xj^y^^z^  axis 

system  of  the  velocity  of  0^^  , the  following  further  relationships  are  easily 
established, 


“l  ° 

cu  - sw  + kj^q 

9 

u 

“ CUj  + SWj^  - 

kq 

V - kp  + hr  ; 

V 

■ ”1  * Vi  • 

"i-'i 

II 

eH 

9 

cw  + su  - hj^q 

9 

w 

* CW^  - SUj^  + 

hq 

When  the  same  system  of  forces  is  reduced  to  a force  through  0 and  a moment 
about  an  axis  through  0 or  to  a force  and  moment  centred  on  0^  , the  follow 
ing  relationships  exist  between  the  components. 

= cX  ■ sZ  ; 

B 

cL  - sN  + kj^Y 

= Y ; 

«i 

B 

M + hZ  - kX 

"l 

• cZ  + sX  ; 

- 

cN  + sL  - hj^Y 

• 

Partial  differentiation  of  these  equations  gives  the  following  conversion 
equations , 


X' 

c^X 

sc(X 

+ z ) 

U 

V 

u 

Z' 

c^Z 

+ 

sc(X 

- z ) 

“l 

u 
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w 

X'  = 

C^x 

+ 

sc(X 

- z ) 

''l 
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w 

w 
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M' 
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- sM 
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+ h(cZ 

1 

M' 

''l 

s 
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w 

+ sM 

u 

+ h(cZ 

1 

M' 
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q 

h(Z 

q 
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w 

Y' 

^1 

= 

Y 

V 

L' 

^1 
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cL 

V 

- sN 
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+ k,Y 

1 V 
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^'l 

s 

cN 

V 

+ sL 

V 

- Vv 
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c^L 

P 

. k^ 

* 

L' 
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c^L 
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N' 

P 

= 

c^N 
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L + k,ik,Y  + c(L  + Y )l  + sisN  - c(L  + N ) - k,  (Y  + N )l 
p 11  Iv  V p/  \r  r p Ir  vf 


,L  ) - h,k,Y  + 
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+ k,Y  - sN 

+ c(L 

1 V 1 1 V 

1 1 V 

1 p p 

1 

,Y  ) - h,k,Y  + 

s-(h,Y 

+ k.L  - sL 

+ c(L 

Ip  1 1 V 

1 1 r 

1 V r 

1 

N + hJh,Y  - c(Y  + N )l  + sisL  + c(L  !■  N ) - h,  (L  + Y )> 
r l(lv  r v/  \p  r p Iv  pi 


These  transformation  relationships  cover  the  effect  of  changing  the  centre  of 
gravity  position  (h  0,  k # 0,  )l  = 0)  as  well  as  the  change  from  one  body-axis 
system  to  another  having  a common  origin,  for  example,  from  aero-body  axes  to  a 
body  axes.  As  such  they  serve  to  illustrate  the  point  already  made  about  the 
complicated  nature  of  the  relationship  between  aero-normal ized  derivatives  and 
coefficient  derivatives  for  a general  body  system  of  axes. 

8 CONCLUDING  REMARKS 

Notwithstanding  the  undoubted  success  that  has  attended  the  use  of  aero- 
dynamic derivative  formulations  of  the  forces  and  moments  in  a variety  of  flight 
dynamic  problems,  it  is  as  well  to  remind  ourselves  that  it  is  incorrect  to 
assume  that  this  is  always  the  case.  For  example,  rapid  movement  of  controls 
(motivators)  as  gust  alleviating  devices  within  an  'active  control'  system  may 
require  the  truly  unsteady  nature  of  the  aerodynamic  phenomenon  to  be  represented 
more  exactly.  Nevertheless,  other  requirements  place  a greater  emphasis  on  the 
development  of  formulations  valid  for  large  disturbance  motion. 

To  meet  this  latter  need  experimental  work  is  in  progress  using  free- 
flight  models  and  on  the  extension  of  the  angle-of-attack  range  of  wind-tunnel 


32 


oscillatory  tests.  In  addition  rotary  rigs  are  being  developed.  In  our  efforts 
in  this  direction  we  should  not  forget  that  the  above  more  fundamental  question 
remains  unanswered,  namely,  under  what  conditions,  specified  quantitatively,  is 
the  customary  derivative  formulation  justified. 
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Table  2 


(1) 

US  symbol* 

%e 

"Zde 

^6e 

Cniy 

Cmq 

CZq 

Si 

"^6r 

^*’6a 

^"6  a 

^"dr 

P 

s 

r 

(2)  ISO  symbol 


Multiplying  factor  to  obtain 
(1)  from  (2)** 


with  ■ c,  ^2  “ ^ *"2  " 


* In  this  table  it  has  been  assumed  that  Cv  and  Cg  are  the 

non-dimensional  forms  of  Xy  and  iy  , but  this  is  not  always 
the  interpretation.  For  example,  Etkin  in  Ref  3,  uses  these  as 
the  equivalent  of  the  second  term  in  the  expressions  for  Xy  and 
Z..  . 


X - 2C^ 

u X 


• &).  ‘ 


- 20^  ^ 
e 


CM 

\3u  / 
N /e 


**  Characteristic  length  Table  1 taken  as  equal  to  for 

longitudinal  derivatives  and  equal  to  I2  ioT  the  lateral 
derivatives . 


Table  3 


(1) 

ISO 

symbol 

X 

Z 

u 

Y 

u 

V 

X 

Z 

w 

w 

X 

Z 

n 

Y 

n 

C 

X 

z 

q 

q 

X. 

z. 

w 

w 

M 

M 

u 

M 

w 

n 

Y 

Y 

P 

Y. 

r 

V 

L 

N 

V 

V 

C 

5 

L 

N, 

5 

5 

M 

q 

M. 

w 

(2)  Old  British  symbol  Multiplying  factor  to 
(R  & M 1801)  obtain  (1)  from  (2) 


In  forming  the  multiplying  factor  it  has  been  assumed  that  the 
normalized  derivatives  in  column  (1)  are  obtained  using  as 
characteristic  length  JIq  “ ^1  longitudinal  and  ilg  ~ ^2 

for  lateral  quantities.  Furthermore  the  practice  of  using  the 
tail-arm  *^0^  always  followed,  so  a check  on  the  defini 

tion  is  advisable  before  using  the  factors. 
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>.  k 1 

‘l*  ‘‘ij 
^0 

i } 


f 

i 

t 

L or  JC 


rate  of  change  of  lift  coefficient  with  angle  of  attack  for  a lifting 
surface,  eg  wing 

value  of  aj^  for  the  fin  of  an  aeroplane 
value  of  a^  for  the  tailplane  of  an  aeroplane 
wing  span 

aerodynamic  mean  chord  of  the  wing 

2 

drag  coefficient  = DApV  S 
lift  coefficient  = L/^pV^S 

2 

rolling  moment  coefficient  = L/^pV  Sb 

2 — 

pitching  moment  coefficient  = M/^pV  Sc 

2 

yawing  moment  coefficient  = N/^pV  Sb 

2 

coefficient  of  component  of  force  along  x-axis  * X/%pV  S 

2 

coefficient  of  component  of  force  along  y-axis  * Y/%pV  S 

2 

coefficient  of  component  of  force  along  z-axis  ■ Z/^^pV  S 
drag  force 

height  above  g;round  of  CG  of  aeroplane 
inclination  of  thrust  line  to  x-axis 

position  of  CG  as  fraction  of  c from  the  leading  edge  of  mean  chord 

displacements  of  the  origins  of  two  systems  of  body  axes  (section  7) 
general  representative  length  used  in  normalization  of  moments  etc. 

particular  representative  lengths 

tail  arm  - distance  from  the  wing-body  aerodynamic  centre  to  tailplane 
aerodynamic  centre 

fin  arm  - distance  from  CG  to  fin  a.c 
tailplane  arm  - distance  CG  to  tailplane  a.c 
rolling  moment  (moment  about  x-axis) 

aero-normalized  derivative  of  rolling  moment  with  respect  to  rate  of 
roll,  p 

aero-normalized  derivative  of  rolling  moment  with  respect  to  rate  of 
yaw,  r 

aero-normalized  derivative  of  rolling  moment  with  respect  to  sideslip 
velocity,  v 

aero-normalized  derivative  of  rolling  moment  with  respect  to  aileron 
angle,  C 

aero-normalized  derivative  of  rolling  moment  with  respect  to  rudder 
angle,  5 


k 
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M 

M 

M 

q 


M 

w 


M. 

w 


N 

N 

P 


N 

r 


N 

V 


N 


C 


N 


P 

q 

r 

Q 

s 

s 

t 

t 

T 


u 

V 

w 

V 

X,  y,  z 
z 
X 
X 

q 


LIST  OF  SYMBOLS  (continued) 

Mach  number 

pitching  moment  (moment  about  y-axis) 

aero-normalized  derivative  of  pitching  moment  with  respect  to  rate  of 
pitch 

aero-normalized  derivative  of  pitching  moment  with  respect  to 
velocity  u 

aero-normalized  derivative  of  pitching  moment  with  respect  to 
velocity  w 

aero-normalized  derivative  of  pitching  moment  with  respect  to  rate  of 
change  of  w 

yawing  moment  (moment  about  z-axis) 

aero-normalized  derivative  of  yawing  moment  with  respect  to  the  rate 
of  roll,  p 

aero-normalized  derivative  of  yawing  moment  with  respect  to  the  rate 
of  yaw,  r 

aero-normalized  derivative  of  yawing  moment  with  respect  to  the  side- 
slip velocity,  v 

aero-normalized  derivative  of  yawing  moment  with  respect  to  aileron 
angle,  C 

aero-normalized  derivative  of  yawing  moment  with  respect  to  rudder 
angle,  c 

rate  of  roll,  angular  velocity  component  about  x-axis 
rate  of  pitch,  angular  velocity  component  about  y-axis 
rate  of  yaw,  angular  velocity  component  about  z-axis 
(with  subscript),  relative  kinetic  pressure 
wing  area  (more  generally  a reference  area) 
fin  area 
tailplane  area 
time 

thrust  force 

component  of  V along  x-axis 

component  of  V along  y-axis 

component  of  V along  z-axis 

resultant  velocity  of  aeroplane 

coordinates  in  the  system  of  axes 

(with  subscript)  fin  arm  with  respect  to  x-axis 

component  of  force  along  x-axis 

aero-normalized  derivative  of  X-force  with  respect  to  the  rate  of 
pitch,  q 
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LIST 

OF 

SYMBOLS 

(concluded) 

i 

1 X 

; u 

aero- norma  1 ized 

derivative 

of 

X-force 

with 

respect 

to 

the 

velocity,  u 

^ X 

w 

aero - norma 1 i zed 

derivative 

of 

X-force 

with 

respect 

to 

the 

velocity,  w 

Y 

component  of  force  along  y 

-axis 

Y 

D 

aero-normal ized 

derivative 

of 

Y-force 

with 

respect 

to 

the 

rate  of  roll,  p 

r 

\ Y 

r 

aero -normal ized 

derivative 

of 

Y-force 

with 

respect 

to 

the 

rate  of  yaw,  r 

Y 

aero -normal ized 

derivative 

of 

Y-force 

with 

respect 

to 

the 

sideslip 

V 

velocity,  v 

Z 

component  of  force  along  z 

-axis 

\ 

Z 

q 

aero -normal ized 

derivative 

of 

Z-force 

with 

respect 

to 

the 

rate  of  pitch. 

Z 

u 

aero-normal ized 

derivative 

of 

Z-force 

with 

respect 

to 

the 

velocity,  u 

1 

■ z 

w 

aero-normal ized 

derivative 

of 

Z-force 

with 

respect 

to 

the 

velocity,  w 

Z 

aero -normal ized 

derivative 

of 

Z-force 

with 

respect 

to 

the 

pitch  motivator 

angle,  n 

a angle  of  attack 

a angle  of  attack  of  a body  x-axis  with  respect  to  equilibrium  flight  path 
e 


Oq  angle  of  attack  of  some  datum  line,  eg  no  lift  line 

8 angle  of  sideslip 

Y angle  of  climb 

r dihedral  angle 

c downwash  angle 

n elevator  angle 

tailplane  setting  angle 
A sweep  angle 

4 aileron  angle 

p air  density 
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i aerodynamic-body  axes 


Fig  4 Attitude  in  pitch  and  angle  of  attack 
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Fig  5 Relationship  between  angle  of  attack  and  attitude  angle 


Fig  6 Normal  velocities  induced 
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Fig  8 Forces  and  moments  due  to  angle  of  attack 
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of  the  associated  moments  (negative) 
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transformations  of  section  7 
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